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Abstract 

The recent design of a nanoscale device with a memristive characteristic has had a 
great impact in nonlinear circuit theory. Such a device, whose existence was predicted 
by Leon Chua in 1971, is governed by a charge-dependent voltage-current relation of 
the form v = M{q)i. In this paper we show that allowing for a fully nonlinear charac- 
teristic v = i](q, i) in memristive devices provides a general framework for modeling and 
analyzing a very broad family of electrical and electronic circuits; Chua's memristors 
are particular instances in which rj(q, i) is linear in i. We examine several dynamical 
features of circuits with fully nonlinear memristors, accommodating not only charge- 
controlled but also flux-controlled ones, with a characteristic of the form i = (,(ip,v). 
Our results apply in particular to Chua's memristive circuits; certain properties of 
these can be seen as a consequence of the special form of the elastance and reluctance 
matrices displayed by Chua's memristors. 
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1 Introduction 



In 1971 Leon Chua predicted the existence of a fourth basic circuit element which would be 
governed by a flux-charge relation, having either a charge-controlled form cp = <f)(q) or a flux- 
controlled one q = a (ip) [10]. This characteristic was somehow lacking in electrical circuit 
theory, since resistors, capacitors and inductors are defined by voltage-current, charge- voltage 
and flux-current relations, respectively. The design of such a memory-resistor or memristor 
at the nanometer scale announced by an HP team in 2008 [37] has driven a lot of attention 
to these devices. Many applications are being developed (cf. [9j [2TJ [361 E3 HQ] and references 
therein). Memristive devices pose challenging problems at the device and circuit modeling 
levels [H [281 H3J [39] but also from a dynamical perspective [201 USE EH (27J [31] . 
The memristor reported in [37] is governed by a relation of the form 

v = M(q)i, (1) 

being framed in the charge-control setting postulated by Chua in his seminal paper [TO] ; 
indeed, assuming to be C 1 , taking time derivatives in the equation ip = <p(q) and using the 
relations <p' — v , q' — i, we are led to (jTJ) with M(q) = (f)'(q). Chua and Kang considered in 
[TT] the more general characteristic 

v = M(q,i)i. (2) 

Both (TjQ) and (T5|), but also the maps governing other devices such as capacitors or resistors, 
are particular instances of the fully nonlinear relation 

v = T](q,i) (3) 

to be discussed in this paper. The general characteristic ()3|) will give a deeper insight into 
the mathematical properties that underly several dynamical features of memristive circuits. 
In particular, memristors increase the dynamical degree of freedom regardless of their ac- 
tual location in the circuit, contrary to capacitors or inductors which in certain (so-called 
topologically degenerate) configurations do not increase the state dimension; additionally, 
memristors have been observed to introduce null eigenvalues in the linearization of differ- 
ent circuits [251 EI]- This kind of problems can be nicely addressed in the framework here 
introduced. 

This will be possible because of the fact that ([3]), together with its dual relation 

i = C(<P,v) (4) 

and their time-varying counterparts, allow for a surprisingly simple, albeit general, model 
for the dynamics of nonlinear circuits including resistors, memristors, capacitors, inductors, 
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and voltage and current sources. Indeed, such a model can be written as 



q' = i g (5a) 

<f' = v v (5b) 

v q = f(q,i q ,t) (5c) 

H = 9(<P,v<p,t) (5d) 

= B q v q + B^v^ (5e) 

= D q i q + D v i v , (5f) 



where B = (B q B v ) and D = (D q D v ) stand for the reduced loop and cutset matrices of the 
circuit (cf. subsection 13 . 1 1) . The circuit elements (and hence B and D) are divided into so- 
called q- and v?-devices, as detailed in Section [2j The key aspect supporting the model ([5]) is 
that e.g. (I5dl) accommodates a great variety of devices, including flux-controlled memristors, 
linear inductors, Josephson junctions, pn and tunnel diodes, independent current sources, 
voltage-controlled current sources, etc.; similar remarks apply to ( l5cl) . 

The goal of our research is two- fold: first, from a modeling point of view, we want not only 
to accommodate in circuit theory general devices governed by ([3]) or (B|, but also emphasize 
the fact that these characteristics allow for a very general description of the circuit dynamics, 
providing in turn a framework for the analysis of different nonlinear phenomena; from an 
analytical standpoint, we wish to tackle several dynamical properties of circuits with fully 
nonlinear memristors. The structure of the paper reflects the two main aspects that drive 
this research: we address modeling issues in Section |2} and dynamical properties in Section 
|3j Section H] compiles some concluding remarks. 

2 Device models 
2.1 g-devices 

Definition 1. A q-device is a circuit element which admits a C l description of the form 

v = r](q,i,t). (6) 

The terminology comes from the fact that g-devices may be controlled by the charge q (this 
will be the case for capacitors), its time derivative q' = i (for current-controlled resistors) 
or both (for q- memristors) . Voltage sources will be included in this group for the sake of 
completeness. 

The relation ([6]) is a time- varying generalization of ([3]). Depending on the actual form of 
the map rj in , g-devices particularize to the ones listed in Table [TJ It is worth noting that 
r) need not be a scalar map; this accommodates coupling effects within the different types 
of g-devices. Keep also in mind that e.g. drj/dq ^ (resp. = 0) means that this derivative 
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does not vanish (resp. does vanish) identically; in the first case, it may find a zero at certain 
points, though; for instance, it vanishes at i = in Chua's memristor, for which v = M(q)i. 



1. A q-memristor is a g-device for which 

The map governing the set of g-memristors will be denoted by v m = r)i(q m ,i m ,t). 

2. A capacitor is a g-device for which 

oq oi 

The relation governing the set of capacitors will be denoted by v c = r)2(q c ,t). 

3. A current- controlled resistor is a g-device for which 

oq oi 

The map governing the set of current-controlled resistors will be written as v r = r] 3 (i r ,t). 

4. An independent voltage source is a g-device for which 

?R = o ?R = o 

dq ' di 

The relation governing independent voltage sources will be denoted by v u = rj^t). 



Table 1: g-devices. 

In a time-invariant setting, if 77 is linear in % then we get Chua's memristor, for which 
v = M(q)i. For this reason we use the expression "fully nonlinear" to label memristors with 
the general form v = rj(q,i). This should cause no misunderstanding with the nonlinear 
nature of the original flux-charge relation (p = 0(g) in Chua's setting. Note also that these 
memristors are referred to in the literature either as charge-controlled memristors or as 
current- controlled memristors; with the expression "g-memristor" we want to emphasize 
that the charge q is the dynamic variable actually involved in the description of the device, 
as it will happen with the flux <p in y9-memristors. 
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2.2 (^-devices 

Definition 2. A tp -device is a circuit element which admits a C 1 description of the form 

i = (((p,v,t). (7) 
The different types of (^-devices are enumerated in Table EJ 



1. A (p-memristor is a <p-device for which 

dip ' dv 

The map governing the set of (^-memristors will be denoted by i w = (\((p w , v w , t). 

2. An inductor is a <y2-device for which 

Off ov 

The relation governing the set of inductors will be denoted by i\ = (2(^1, t). 

3. A voltage- controlled resistor is a v?-device for which 

Oip ov 

The map governing the set of voltage-controlled resistors will be written as i g = £3(1^, t). 

4. An independent current source is a y?-device for which 

The relation governing independent current sources will be denoted by ij = (^(t). 



Table 2: ^-devices. 

Again, Chua's flux-controlled memristor [TU], for which i = W((p)v, is obtained in par- 
ticular when ( is time-independent and linear in v. 

We do not impose any condition on the time derivatives, either for q- or ^-devices. If they 
vanish in memristors, resistors, capacitors and inductors we would be led to a time-invariant 
setting, with r] k , ( k {k = 1, 2, 3) being independent of t. For sources this assumption would 
model DC ones. It is also worth remarking that current-controlled voltage sources (CCVS's) 



5 



could be easily included in the above taxonomy as g-devices and, similarly, voltage-controlled 
current sources (VCCS's) can be modeled as y?-devices. 



2.3 Characteristic matrices and passivity 

The matrices of partial derivatives of the maps introduced above define the so-called char- 
acteristic matrices of the different devices. 

Definition 3. The incremental memristance, elastance and resistance matrices of q-memristors, 
capacitors and current- controlled resistors, respectively, are defined as 

M(q m , % m , t) = — , E c (q c , t) = — — , R(i r , t) = — — . 

Additionally, the incremental elastance of q-memristors is 

dq m 

In Chua's memristor, for which v = M(q)i, the incremental memristance depends only 
on the charge q. The fact that q(t) = i m (r)dr makes this element behave as a resistor 
in which the resistance depends on the device history; the "memory-resistor" name stems 
from this. 

If the elastance of capacitors is non-singular at a given operating point, then because 
of the implicit function theorem these devices admit (at least locally) a voltage-controlled 
description 

q c = 7(u c ,f). 

The incremental capacitance matrix is then defined as 

C(v c ,t) = = (EMvc,t),t))-\ 

Note that when a given set of devices (memristors, capacitors or resistors) does not exhibit 
coupling effects, the corresponding characteristic matrix is a diagonal one, with diagonal 
entries defined by the memristances, elastances or resistances of the individual devices. 

The characteristic matrices of y?-devices are defined analogously, as detailed below. 

Definition 4. The incremental memductance, reluctance and conductance matrices of tp- 
memristors, inductors and voltage-controlled resistors, respectively, are defined as 

ur( +\ d (i(¥w,v w ,t) dC 2 ((pi,t) d( 3 {v g ,t) 

W(ip w ,v w ,t) = , Tli{(pi,t) = — , G(vg,t) = — — • 

ov w dipi dv g 

The incremental reluctance of ip-memristors is 

v , ,N dCl(<Pw,Vw,t) 

rC w yP w ,V w ,T) 
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As for capacitors, if the reluctance of inductors is non-singular, then via the implicit 
function theorem these devices can be shown to admit a local current-controlled description 

<Pi = £(ii,t), 

and the incremental inductance matrix is defined as 

In order to define the local notions of passivity and strict passivity, we make use of the 
concept of a positive (semi) definite matrix; a square matrix P is positive semidefinite (resp. 
definite) if u T Pu > (resp. > 0) for every non-vanishing vector u. We do not assume P to 
be symmetric. 

Definition 5. We call the q-memristors, capacitors, current- controlled resistors, ip-memris- 
tors, inductors or voltage-controlled resistors locally passive (resp. strictly locally passive,) at 
a given point if the incremental memristance, elastance, resistance, memductance, inductance 
or conductance matrix is positive semidefinite (resp. definite) at that point. 

It is easy to check that a positive definite matrix is non-singular, and that its inverse is 
itself positive definite; this means that the capacitance and inductance matrices are positive 
definite if the elastance and reluctance matrices of capacitors and inductors are so. No 
condition is imposed on the elastance and reluctance matrices of q- and <y2-memristors, except 
for the fact that they cannot vanish identically since otherwise these devices would amount 
to (current- or voltage-controlled) resistors. Note finally that, with terminological abuse, the 
adverb "locally" is often omitted when a given set of devices satisfies the (semi)definiteness 
requirement everywhere. Similarly, the term "incremental" is very often omitted in the 
characteristic matrices. 

3 Analytical results: nondegeneracy, null eigenvalues 

The general framework introduced in Section [2] makes it possible to address different analyti- 
cal properties of memristive circuits in broad generality. After introducing some background 
material in subsection I3.1[ we tackle in subsections 13. 2\ 13.31 and 13.41 two particular problems 
which rely on linearization, emphasizing the role of the form of the different characteristic 
matrices discussed above; specifically, the absence of memristance (resp. memductance) in 
capacitors (resp. inductors), namely, the conditions drj/di = 0, d(/dv = holding for them, 
together with the form of the elastance and reluctance matrices for different memristors, 
will explain certain dynamical features of memristive circuits. Our approach should also be 
useful in future analyses of many other aspects of memristive circuit dynamics, including 
e.g. stability properties, oscillatory phenomena or bifurcations. 
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3.1 Some auxiliary results from digraph theory 



Many aspects of our analysis will rely on the properties of the directed graph or digraph 
underlying a given electrical circuit. We compile below some background material on digraph 
theory. The reader is referred to [U |2j [5j [HI [30] for additional details. 

We will work with a directed graph having n nodes, m branches and k connected compo- 
nents. A subset K of the set of branches of a digraph is a cutset if the removal of K increases 
the number of connected components of the digraph, and it is minimal with respect to this 
property, that is, the removal of any proper subset of K does not increase the number of 
components. In a connected digraph, a cutset is just a minimal disconnecting set of branches. 
The removal of the branches of a cutset increases the number of connected components by 
exactly one. Furthermore, all the branches of a cutset may be shown to connect the same 
pair of connected components of the digraph which results from the deletion of the cutset. 
This makes it possible to define the orientation of a cutset, say from one of these components 
towards the other. The cutset matrix D = (d^) is then defined by 



The rank of D can be proved to be n — k; any set of n — k linearly independent rows of D 
defines a reduced cutset matrix D G ^(n-k)xm^ i n a connected digraph, any reduced cutset 
matrix has order [n — 1) x m. 

The space spanned by the rows of D equals the one spanned by the rows of the so-called 
incidence matrix, defined as A = (a^-) with 



Similarly, chosen an orientation in every loop, the loop matrix B is defined as (bij), where 



The rank of this matrix equals m — n + k. A reduced loop matrix B is any ((m — n + k) x m)- 
submatrix of B with full row rank. 

We denote by B^ (resp. Bg_x) the submatrix of B defined by the columns which cor- 
respond to branches belonging (resp. not belonging) to a given set of branches K; the same 
applies to the cutset matrix D. Certain submatrices of B and D characterize the existence of 
so-called f^-cutsets and i^-loops (that is, cutsets or loops just defined by branches belonging 
to K), as stated below (cf. [30l Sect. 5.1]). 




1 if branch j is in cutset i with the same orientation 
1 if branch j is in cutset i with the opposite orientation 
if branch j is not in cutset i. 




1 if branch j leaves node i 
1 if branch j enters node i 
if branch j is not incident with node i. 




1 if branch j is in loop i with the same orientation 
1 if branch j is in loop i with the opposite orientation 
if branch j is not in loop i. 



S 



Lemma 1. Let K be a subset of branches of a given digraph Q . The set K does not contain 
cutsets if and only if Bk has full column rank or, equivalently, iff Dq-k has full row rank. 

Analogously, K does not contain loops if and only if Dk has full column rank or, equiv- 
alently, iff Bg_ K has full row rank. 

Actually, the dimensions of the spaces ker Bx and ker are defined by the number of 
linearly independent K-cutsets and if-loops, respectively. 

The proof of the following result can be found e.g. in [T4"l Sect. 7.4]. 

Lemma 2. If the columns of the reduced loop and cutset matrices B, D of a digraph are 
arranged according to the same order of branches, then BD T = 0, DB T = 0. 

The relations imD T = kerB and imB T = kerD hold true. This expresses that the cut 
space imD T spanned by the rows of D can be described as ker B and, analogously, the cycle 
space imB T spanned by the rows of B equals kerD [5]. These spaces are orthogonal to each 
other since (imD 7 ) 1 - = (kerB) 1 - = imB T . 

Finally, when applying these results to circuit analysis, we will split the columns of B and 
D according to the nature of the devices accommodated in the corresponding branches; B q , 
B v , D q , Dp describe the sub matrices of B and D defined by q- and y?-devices; these matrices 
will be further split into B m , B c , B r , B u , B w , Bi, B g , Bj (resp. D m , D c etc.), according to 
the division of q- and (^-devices into g-memristors, capacitors, current-controlled resistors 
and voltage sources, and (^-memristors, inductors, voltage-controlled resistors and current 
sources, respectively. With this notation, the identity BD T = in Lemma [2] can be written 
as 

B m D T m + B C D T C + B r Dj + B u Dl + B W D T W + B&J + B g Dj + BjDj = 0. (8) 
3.2 Nondegenerate problems and the order of complexity 

The order of complexity or state dimension of an electrical circuit is the number of variables 
that can be freely assigned an initial value. Stemming from the work of Bashkow and Bryant 
[3l [7J [8] , a circuit composed of capacitors, inductors, resistors and voltage and current sources 
is said to be nondegenerate if its order of complexity equals the number of reactive elements 
(capacitors and inductors). A necessary condition for a circuit to be nondegenerate is its 
topological nondegeneracy, that is, the absence of VC-loops (loops formed only by voltage 
sources and/or capacitors) and IL-cutsets (cutsets defined only by current sources and/or 
inductors). Different sufficient conditions can be given, involving e.g. the strict passivity of 
the circuit matrices or the structure of the circuit spanning trees [121 ESS E2J [3H [35] . 

In Theorem [1] we address the characterization of the order of complexity of circuits with 
fully nonlinear memristors, under strict passivity assumptions and restricting the discus- 
sion to cases without VC-loops and IL-cutsets. We show that memristors, even under a 
fully nonlinear assumption, do not introduce topological degeneracies; this means that ev- 
ery memristor increases by one the order of complexity, regardless of its location in the 
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circuit. The key difference with capacitors and inductors, which actually introduce degen- 
eracies when entering VC-loops or IL-cutsets, is made by the fact that the characteristics of 
these, namely, v c = r)2(q c ,t) and %i = (2(^1, t), do not involve the current i c or the voltage vi, 
respectively 

Note, in this regard, that the equations q' r = i r , q' u = i u , (p' g = v g , <p'a = vj will be excluded 
from the dynamical description of the circuit, since the variables q r , q u , ip g , ipj are decoupled 
from the rest of the system and are hence irrelevant from the dynamical point of view. This 
means that the circuit dynamics will be defined by the model 







(9a) 


<Pwl 


= v wt 


(9b) 





Vq /(<7mc? iqi t) 


(9c) 





= i<p- g(<Pwi,v v ,t) 


(9d) 





= B q Vq + ByVy 


(9e) 







(9f) 



in the understanding that q mc and ip w i stand for (q m , q c ) and ((p w , <pi), respectively. The 
same splitting applies to i mc and v w i. Notice that / and g group together the maps r] k , ( k 
(k = 1, ... ,4) arising in Tables [T] and [2] for the sets of q- and (^-devices. 

Theorem 1. Consider a circuit without VC-loops and IL-cutsets in which the memristance, 
resistance, memductance and conductance matrices M , R, W , G are positive definite. Then 
its order of complexity is given by the total number of q- and ip-memristors, capacitors, and 
inductors. 



Proof. We will check that the relations (}9c^-(l9fl) make it possible to write explicitly all the 
branch voltages and currents in terms of q m , q c , <p w , (pi, and that these relations impose 
no constraint among these variables in the absence of VC-loops and IL-cutsets. The way 
to do so is to check that the matrix of partial derivatives of (|9cj)-(l9fl) with respect to the 
branch voltages and currents v q , i q , v 9 , i v , is non-singular; a straightforward application of 
the implicit function theorem yields the result. 
This matrix of partial derivatives has the form 



J 



/ 


h 


— M q 








\ 










-w v 


h 






Bq 












V 





Dq 





D 9 


J 



(10) 



with 



10 



Using a Schur reduction [191 EO], it is easy to see that ( TlQjl is non-singular if and only if so 
it is 



According to the classification of g-devices into g-memristors, capacitors, resistors and 
voltage sources, the matrix M q has the block-diagonal structure block-diag{M, C , R, U }, 
where M and R are the incremental memristance and resistance matrices. Analogously, 
reads as block-diagjW 7 , 0;, G, 0j}. This confers the matrix J red in f lTTj) the form 



B m M B r R B w B[ Bg Bj 
D m D c D r D u D W W D g G 



;i2) 



Assume that ( 1T2|) has a non-trivial left-kernel; that is, suppose that there exists a non- 
vanishing (x T y T ) such that the following relations hold: 



x T B m M + y T D m 


= 


(13a) 


y T D c 


= 


(13b) 


x T B r R + y T D r 


= 


(13c) 


y T D u 


= 


(13d) 


x T B w + y T D w W 


= 


(13e) 


x T B l 


= 


(13f) 


X T Bg + y T DgG 


= 


(13g) 


x T B ] 


= 0. 


(13h) 



Multiply the identity OH]) from the left by x T and from the right by y. Using f ]13bj) . (113dj) . 
f fT3f|) and f fT3hj) . we get 

x T {B m Dl + B r D T r + B W D T W + BgD T g )y = 0. 

This equation can be recast, using (I13ap . (I13cj) . (I13ep and (|13gP, as 



x T B m M T B T m x + x T B r R T Bjx + y T D w WD T w y + y T D g GD T g y = 0. (14) 

Be aware of the fact that M T and R T are positive definite since so they are M and R. 
Together with the positive definiteness of W and G, this implies that the relations 

x T B m = (15a) 

x T B r = (15b) 

y T D w = (15c) 

y T D g = (15d) 
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follow from (1141) . These relations yield, in turn, 



y T D m = (16a) 

y T D r = (16b) 

x T B w = (16c) 

x T B g = 0. (16d) 

Now, the identities (Jl3f|) . f ll3hl) . ffTHal) . f lTobl) . ffTocl) and f ITMl) imply that x = 0, because 
of the absence of VC-loops and Lemma [U Similarly, ( 113bl) . ( 113dl) . (I15cf) . ( 115dl) . (I16a.j) and 



( 116bl) . together with Lemma [Hand the absence of IL-cutsets, yield y = 0. This means that 
( TP21) and so J in ( TTOl) are indeed non-singular; therefore, the implicit function theorem makes 
it possible to write v q , v v , i q , i v in terms of q m , q c , (p w and <pi, and this in turn yields an 
explicit ODE formulated in terms of the charges of g-memristors and capacitors and the 
fluxes of yj-memristors and inductors. The state dimension of this ODE is obviously defined 
by the number of memristors, capacitors and inductors and the proof is complete. □ 



Note that this result applies to fully nonlinear memristors, regardless of the actual form 
of the maps r/i and (1, as far as the memristance M and the memductance W are positive 
definite. It holds true, in particular, for Chua's memristors. 

The different role played in this regard by g-memristors and capacitors, on the one hand, 
and by (^-memristors and inductors, on the other, becomes apparent in the light of the form 
of the matrix f lT2|) . Indeed, the fact that the characteristic v c = r] 2 {q C it) of capacitors does 
not involve the current i c is responsible for the vanishing block over D c in (|T2|) . Together 
with the null block over D u , this explains the key role of VC-loops in the state dimension 
problem. By contrast, the presence of the M matrix within ( !T2|) makes the location of mem- 
ristors irrelevant in this problem, showing that g-memristors cannot introduce topological 
degeneracies. In particular, VCM-loops (with at least one g-memristor) do not reduce the 
state dimension of the problem. The same reasoning applies to yj-memristors, inductors, 
IL-cutsets and ILW-cutsets. 

With more technical difficulties it is possible to show, assuming the elastance of capacitors 
and the reluctance of inductors (or, equivalently, the capacitance and the inductance) to 
be positive definite, that the order of complexity of topologically degenerate circuits with 
fully nonlinear memristors is given by the total number of memristors plus the number 
of capacitors in a normal tree (that is, a tree including all voltage sources, the maximum 
possible number of capacitors, the minimum possible number of inductors and no current 
source) and the number of inductors in a normal cotree. This provides an extension of 
Bryant's results [7J [8] to circuits with fully nonlinear memristors. 
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3.3 Null eigenvalues; regular equilibrium points 



Many qualitative properties of dynamical systems at equilibria can be characterized in terms 
of the linearized problem. In particular, zero eigenvalues of the linearization may be respon- 
sible for stability changes and bifurcation phenomena. In particular, within the context of 
memristive systems, some particular circuits with Chua-type memristors have been shown 
in the literature to display null eigenvalues [25] EE] • Below we tackle this problem for general 
circuits with fully nonlinear memristors, addressing Chua's memristors as a particular case 
with distinctive properties. 

System is a semiexplicit differential-algebraic equation (DAE) [6] [22] [23], [291 130] . 
Equilibrium points are defined by the vanishing of its right-hand side, and the linearization 
at a given equilibrium leads to the matrix pencil [[15] XH — K, where H = block-diag{/, 0} 
and K is the Jacobian matrix of the right-hand side at equilibrium. The spectrum of the 
matrix pencil is the set of values of A which make XH — K singular. In particular, the pencil 
has null eigenvalues if and only if the Jacobian matrix K is singular. 

An equilibrium point is said to be regular if and only if if is a non-singular matrix. 
Regular equilibria are important regarding DC-solvability and Newton-based computations 
[IB"] . The following result extends to systems with fully nonlinear memristors a property 
already known for RLC circuits [16], [17] [24] . With the same terminological convention, 
a VLW-loop is a loop defined by voltage sources, inductors and/or yj-memristors, and an 
ICM-cutset is a cutset including only current sources, capacitors and/or g-memristors. 

Theorem 2. Consider a memristive circuit with positive definite incremental resistance and 
conductance matrices R, G. An equilibrium point of this circuit is regular if and only if 

• the elastances E m , E c of q-memristors and capacitors, as well as the reluctances TZ W , 
TZi of if -memristors and inductors, are non-singular at the equilibrium; and 

• the circuit does not have either VLW-loops or ICM-cutsets. 
Proof. The Jacobian matrix K reads as 



with 
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and 
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It is clear that the elastances E m , E c and the reluctances 1Z W , TZi must be non-singular 
for K to be so; if these requirements are met, the non-singularity of K relies on that of 
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and, by means of a Schur reduction, on the non- singularity of 

( B m B c B r R B g B 3 
red V D r D u D g G D w D l 

The blocks comprising B m , B c , Bj show, in the light of Lemma (TJ that the absence of ICM- 
cutsets is necessary for K re d to be non-singular; analogously, the blocks including D u , D w , 
Di rule out the presence of VLW-loops. Conversely, in order to show that the absence of 
these configurations guarantees K re & to be non-singular, it suffices to check that a column- 
reordering gives this matrix the form of J re d in ( |TT|) ; proceeding exactly as in the proof of 
Theorem [T], and using only the positive definiteness of R and G, it is easy to check that 
in the absence of VLW-loops and ICM-cutsets the matrix K re d (and then K) is actually 
non-singular, as we aimed to show. □ 
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More is true, actually; proceeding as in Theorem [3] below, one can show that, in problems 
in which the elastance of g-memristors and capacitors and the reluctance of 99-memristors 
and inductors are non-singular, the geometric multiplicity of the zero eigenvalue of K equals 
the number of independent VLW-loops and independent ICM-cutsets, where the notion of 
independence relies on the cut- and cycle-spaces described in subsection 13.11 

3.4 Null eigenvalues in circuits with Chua-type memristors 

A distinct feature of Chua's memristors is that the elastance E m = dr]i/dq m and the re- 
luctance TZ W = d(i/dtp w do vanish at equilibrium points, because T]i(q m , i m ) = M(q m )i m 
and (i((p w , v w ) = W((p w )v w are linear in i m and v w , respectively, and the identities i m = 0, 
v w = hold at equilibria. Theorem |2] above then predicts the existence of null eigenvalues 
in Chua-type memristive circuits owing to this zero-crossing property, as already observed 
in specific examples [25l l3T]. The number of memristive devices actually characterizes the 
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geometric multiplicity of the zero eigenvalue in a broad class of circuits with Chua-type 
memristors, as detailed below. 

Theorem 3. Consider a circuit in which the elastance E m of q-memristors and the re- 
luctance 1Z W of ip -memristors do vanish at equilibria. If the resistance and conductance 
matrices R, G are positive definite at a given equilibrium, and the elastance E c of capaci- 
tors and the reluctance 1Z[ of inductors (equivalently, the capacitance and the inductance) 
are non-singular, then the geometric multiplicity of the null eigenvalue equals the number of 
memristors, independent VL-loops and independent IC-cutsets. 

Proof. Making E m = and 1Z W = in ( fl8l) . the corank of the matrix K can be easily seen 
to equal the number of (q- and (p-) memristors plus the corank of 
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(20) 



The corank of (1201) in turn equals that of 

/ B C E C B r R B g Bj \ , . 

V D{Ri D r D u D g G / ' 1 ' 

Contrary to (fl9l) . ( J2TT) is not a square matrix, and its corank is defined by the dimension of 
its (right) kernel. A vector (u,v,w,x,y, z) of the kernel of this matrix verifies 
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The relations ([22]) imply that (0, E c u, Rw, 0, 0, 0, y, z) e ker (B m B c B r B u B w B { B g Bj) = 
keiB, and (0, 0, w, x, 0, IZiv , Gy, 0) G ker(D m D c D r D u D w D\ D g Dj) = kerD. Use then 
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the identities ker B = imD T and kerD = imB T to derive the existence of vectors p, q 
satisfying 

= Dip, E c u = D T c p, Rw = D T r p, = D T u p, = D T w p, = D?p, y = D T g p, z = Djp, 
= Big, = B T c q, w = Bjq, x = B T u q, = B%q, 7L x v = B?q, Gy = B T g q, = Bjq. 

Pre- and post-multiplying (jSj) by q T and p, respectively, we get, in the light of these identities, 

w T Rw + y T G T y = 0, 

that is, w = y = because of the positive defmiteness of R and G. Then (E c u, z) belongs to 
ker(S c Bj) and (x,1Ziv) belongs to ker(D u D{). Due to the non-singular nature of E c and 
Hi, this means that the corank of (120]) equals the number of independent IC-cutsets plus the 
number of independent VL- loops (cf. subsection 13. lft . Hence, the corank of the matrix K, 
and therefore the geometric multiplicity of its null eigenvalue, equals the number of (q- and 
if-) memristors plus the number of independent VL-loops and the number of independent 
IC-cutsets, as we aimed to show. □ 

In particular, in circuits without VL-loops and IC-cutsets, the geometric multiplicity 
of the null eigenvalue matches exactly the number of memristors; this is the case in the 
examples arising in [251 13T] . 

4 Concluding remarks 

The fully nonlinear characteristics (EJ) and (jlj) introduced in the present paper might accom- 
modate future devices arising in nonlinear circuit theory and displaying memristive effects, 
or provide more accurate models for already existing devices. The general form of these 
relations and their time-varying counterparts allows for a simple description of a very broad 
class of electrical and electronic circuits; different dynamical features of these can be then 
addressed in great generality. Several of these features can be understood to follow from 
the special form of the characteristic matrices of the different devices arising in the analy- 
sis. For instance, the zero-crossing property of Chua-type memristors makes them display a 
vanishing elastance and reluctance at equilibria; from Theorems Q] and [3j it then follows that 
in strictly passive circuits with Chua-type memristors, not exhibiting VC-loops, IL-cutsets, 
VL-loops or IC-cutsets, the order of complexity equals the number of memristors plus reac- 
tive elements, every memristor introducing a vanishing natural frequency in the linearized 
problem. Similarly, the absence of null eigenvalues is characterized, for general memristive 
circuits and in terms of the circuit matrices and the digraph topology, in Theorem [2j The 
characterization of other analytical properties of memristive circuits within the framework 
here introduced defines the scope of future research. 
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